Stochastic resonance in a suspension of magnetic dipoles under shear flow. 
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, We show that a magnetic dipole in a shear flow under the action of an oscillating magnetic field 

' displays stochastic resonance in the linear response regime. To this end, we compute the classical 

£\j | quantifiers of stochastic resonance, i.e. the signal to noise ratio, the escape time distribution, and 

■ the mean first passage time. We also discuss limitations and role of the linear response theory in its 

' applications to the theory of stochastic resonance. 
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I. INTRODUCTION 

1 ; The < _ b ee„ . _ M < _ _ . ,, t _ g. 

■ This kind of systems arise commonly in fields of physics, chemistry and biology. Examples are found in problems 
involving transport at the cellular level optical and electronic devices j| and signal transduction in neuronal 

tissue [||J|], to cite just a few. 

A particularly interesting phenomenon, occurring in periodically-driven non-linear noisy systems, is stochastic 
resonance (SR) This term refers to the enhancement of the response of the system to a coherent signal when 
the intensity of the noise grows, instead of the degradation that one naively expects. The mechanism leading to this 
phenomenon is quite simple. Imagine a system which exhibits an energetic activation barrier. In the presence of noise, 
the system can be assumed to surmount this barrier with a rate proportional to e~~ AE / D , where AE is the height 
I ' of the barrier and D is the intensity of the noise acting on the system. The inverse of this rate defines the average 
waiting time, T(D), between two noise-induced transitions. In the presence of a periodic forcing, the height of the 
barrier is periodically raised an lowered. When the period of the external force synchronizes with 2T(D), the barrier 
t***" \ surmounting will be enhanced by the cooperative effect of the noise and the periodic forcing. 

Although originally proposed for systems in a double- well potential, this original scheme has been overcome. In fact, 
it is known that SR is exhibited by several classes of monostable systems, among which one could mention excitable 
and threshold systems ]^[Tl|l or systems which do not follow an activated dynamics but a relaxational dynamics 

us 

In this paper we will show that a magnetic dipole immerse in a shear flow exhibits stochastic resonance when a 
weak oscillating magnetic field is acting on it. The presence of this flow takes the system out of equilibrium impressing 
certain peculiarities in the behaviour of this system. In order to treat this problem we will analyze the response of 
the system in the linear regime. Previous study of the dynamics of a dipole under an oscillating magnetic field has 
revealed that linear response theory (LRT) predicts a monotonically decreasing behaviour for the ratio between the 
output signal and the output noise or signal to noise ratio (SNR), i.e. for very weak applied fields the dipole does 
not exhibit SR [fijl. In the present case there is a new ingredient, absent in fjlJ]: the presence of the shear flow, 
which is determinant for many interesting aspects of the dynamics of this system. Additionally, although we show 
the existence of SR in the linear regime, we discuss limitations and role of the LRT in its application to the theory of 
SR, mainly related to questions concerning the fluctuation-dissipation theorem. 

The paper is organized as follows: in section II we analyze the dynamics of a dipole in a shear flow and find the fixed 
points. Section III is devoted to study the response of the system to an oscillating magnetic field by computing the 
susceptibility. In section IV we calculate the power spectrum and the signal-to-noise ratio. In section V we compute 
the escape time distribution and from it the mean first passage time. Finally, in section VI we discuss our results. 



II. DYNAMICS OF A DIPOLE IN A SHEAR FLOW: FIXED POINTS AND THEIR STABILITY 

We consider a dilute colloidal suspension of ferromagnetic dipolar spherical particles, with magnetic moment, 
rfi = m s R, where R is an unit vector accounting for the orientation of the dipole, therefore, being the magnetic 
moment rigidly attached to the particles. Each dipole is under the influence of a shear flow with vorticity Q — 2ujqz, 
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with z being the unit vector along the z-axis, and of an oscillating field H = He iut x, with x being the unit vector 
along the x-axis. The dynamics of these dipoles is governed by the following equation of motion 

where / is the moment of inertia of the particles, £ r — 8iri]Qa 3 is the rotational friction coefficient, with 770 the solvent 
viscosity, and a the radius of the particle. For t 3> 7y,with r r = I/£,r being the inertial time scale, the motion of 
the particle enters in the overdamped regime. This time scale defines a cut-off frequency uj r — r~ , such that the 
condition for overdamped motion is equivalent to w < oj r . In this case eq. ([!]) yields the balance condition between 
the magnetic and hydrodynamic torques acting on each particle 

mx H + £ r (J^n-ttp^j = 0, (2) 

which, together with the rigid rotor evolution equation 



dR 



— =QpX R, (3) 



leads to the dynamic equation for R 



dR 
~dt 



u {t + A(f) (R x X x R. (4) 



Here A(i) = (m s H / ^ r ujo)e~' luJt , with fl p being the angular velocity of the particle. 

The computation of the fixed points of eq. (Q) when the magnetic field is held constant, i.e. X(t) = Xq = cnt., and 
their linear stability analysis are given in detail in Appendix A. After some algebra eq. ([|) becomes 



dR 

-dt= U ° 



zxR+Xx- XR{R ■ x) . (5) 



For Ao > 1, this equation has only a linearly stable stationary state. The orientation of the suspended particles is 
fixed to 



R s = ^l-X- 2 £ +X~ 1 ~y. (6) 

This means that in this regime the hydrodynamic torque, which tends to cause the rotation of the particles is 
insufficient to overcome the magnetic torque which maintains constant their orientations. 

For Ao < 1, which is the case we arc interested in, the particles undergo a rotation around a fixed axis lying in the 
y-z plane, being the director of this axis given by 



R s =± VI - X 2 ij + Xk (7) 

In this case the hydrodynamic torque is strong enough to make the dipole to precess around the orientation R s eq. 
(0) (see Appendix A). 

III. RESPONSE TO AN OSCILLATING MAGNETIC FIELD 

The analysis of Section II has been carried out for the deterministic dynamics of a magnetic dipole in a shear flow. 
Fluctuations are introduced by means of a Brownian torque. The corresponding Langevin equation is 

(lh cd \x{t)(Rxx)+f+-J— (RxF B (t))\ xR, (8) 



dt [ V / 

where Fg (t) is a Gaussian white noise of zero mean and correlation function 
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(F B (t)F B (t')) = 2i r K B T5(t - t'). (9) 
The Fokker-Planck equation corresponding to eq. (|^) is given by 

dM&,t) = (Co + XWdM^t), (10) 
where Cq and C\ are operators defined by 

C = -uo~z-TZ + D r TZ 2 , 

Ci = 2lj R ■ ~x- uj {R x f) • TZ, (11) 

with D r = k B T/^ r being the rotational diffusion coefficient, and TZ = R x d/dR the rotational operator. Notice 
that the first and second terms on the right hand side of eq. (p"l|)i, correspond to convective and diffusive term, 
respectively. Moreover eq. ( |l0| ) which, according to eq. (0), rules the Brownian dynamics in the case of overdamped 
motion, is valid in the diffusion regime. This regime is also characterized by the condition f > r r , or equivalently 
lu <C Lj r , which implicitly involves the white noise assumption. 

To solve the Fokker-Planck equation ( [l0| ) we will assume that Ao = |A(i)| constitutes a small parameter such that 
this equation can be solved perturbatively. Thus up to first order in A, the solution of the Fokker-Planck equation 
©is 



V(R,t) = e (t -* o)£o * (t )+ / dt'A(i')e (t -*' )£o £i*o(t') (12) 
Here ^o(t') — e*-* — *o)^o nj>o ^ _ ^ - 1S zeT0 order solution at time t' , and 

t = t ) = S(R - R ), (13) 



with Ro being an arbitrary initial orientation. As follows from eq. (|1 1|) i , the unperturbed operator Co is composed of 
the operators TZ Z and TZ 2 , which are proportional to the orbital angular momentum operators of quantum mechanics 
L z and L 2 , respectively, and, therefore, their eigenfunctions are the spherical harmonics fl4| 

Tl z Y lm (R) = imY lm (k), 

TZ 2 Y lm (R) = -l(l + l)Y lm (R). (14) 

Given that we know how TZ acts on the spherical harmonics, it is convenient to expand the initial condition in 
series of these functions, since the spherical harmonics constitute a complete set of functions which are a basis in the 
Hilbert space of the integrable functions over the unit sphere ]l5fl 

*°(ii,to) = S(R-Ro) = ]T ]T Y l * m (Ro)Yi m (R) . (15) 

1=0 m=-l 

Using this expansion in eq. (|l2|), for the first order correction to the probability density, Alt' = — ^>o, we obtain 

oo j ft 

A*(R,t) = J2 E dt'Xit^Y^i^e^^Cre^'^^Y^iR). (16) 

;=0 m=-l , ' to 

Notice that the integral of A^f(R, t) over the entire solid angle is zero, in agreement with the fact that the unperturbed 

solution ^o{R,t) is normalized. 

Since, we are interested in the asymptotic behavior we will set to — * — oo. In this limit, eq. ( [12] ) becomes 

V(h,t) = / dt' X(t')e {t - t,)Co 2 k-$X, (17) 

I J — oo 

where now 
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dt' \(t')e {t - t>)Co 2R-'x, 



and 



-in 



(18) 



(19) 



is the uniform distribution function in the unit sphere. 

From eq. ( |l8| ) the contribution of the AC field to the mean value of the orientation vector R can be obtained as 



R(t) = I dRRAV = — I dt'\{t') I dRRe {t ' t,)Co 2R- . 



(20) 



This equation can be written in the more compact form 



dt'X(t') Xi (t-t'), 



Ri{t) = 

where the response function |l6| has been defined as 

Xi(r) = ^ / dkR t e TCo 2k-^ 
for r > 0. 

By causality, we can write t — > oo in the upper limit of the integral in eq. (|2l|); hence, this equation becomes 

Ri(t) = Xi(w)A(t). 

where Xi( w ) is the generalized susceptibility, which is the Fourier transform of Xi( T ) 



(21) 



(22) 



(23) 



Xt(<*>) 



1 

4tt 



dTe lUT I dRR,e TCa 2R-. 



(24) 



From this equation we obtain the components of the susceptibility 

Xx(u) = | { 



X»M = | { 



2D r 


. (w -«0 




iD 2 + (ui-ui a ) 2 


2D r | 




4D^ + (w+w ) 2 ^ 


L iD 2 + (u>+w ) 2 \ 


(wo— 0)) 




iD'i + iuj-ua) 2 


^ 1 iD'i+iuj-ujoy 2 




■ 2D r 


4Dj? + (w+w ) 2 










+ 



(25) 

(26) 
(27) 

The quantities Xx, and Xyt have poles at u> — ±u>o ± 2D r i. The inverse of the imaginary part of these poles, (2D r )~ 1 , 
defines the Brownian relaxation time. 

IV. POWER SPECTRUM 

In order to discern whether or not SR is present in the relaxation process of the system under consideration we 
compute the power spectrum, which, following the Wiener-Khinchine theorem, is given by the Fourier transform of 

the correlation function Jl7|,|T[] . Since we will take as output signal the projection of R parallel to the magnetic field, 
i.e. R x , we only compute the correlation function of this quantity. 
The correlation function of R x is defined by 
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(R x (t)R x (t + T)\R (t )) = dv duu x v x ^(v,t;u,t + T\R ,t ), 



(28) 



where the initial condition is taken as ^f(R, to) = S(R — Ro). The above quantity can be calculated from the solution 
Fokker-Planck equation simply by recalling the following properties of a Markov process |l7|l . 



*(i;i,ti; ...;v n ,t n ) = h\v 2 , h; ...;v n ,t n ), 

*(i/i,ti|u2,i2; ...;v n ,t n ) = ^(vi,ti\v 2 ,t 2 )- 

where t\ > t% > ... > t n . By combination of these two properties eq. (^8j) becomes 

(R x (t)R x (t + T)\R (t )) = / dito x *0t,t\ito,to) / d€tu x *(€t,t + r|€f,«) 



(29) 



(30) 



To proceed further, we compute in first place the integral over u in eq. (BOJ). From eq. (fL2|), if r > 



t+T 



*(tZ, i + r|u, t) = e T£ °<5(u -v)+ ds \{s)e (t+T ' s)Ca de sCa 5(u - v) 



= $ (u,t + T\v, t) + A^{u,t + T\v, t), 

thus we have 

y duu x ^(u,t + r\v,t) — J duu x 1 i'o(u,t + r\v,t) + J duu x A^(u,t + r\v,i). 
The result of these integrals are (for the detailed derivation, see Appendix B) 

2^ r 



(31) 



(32) 



duuMn,t + r\v,t) =-J -f [e-^-^YiiW + e-^+^Yi-^if) 
duu x A^(u, t + t\v, t) = 

/=0 in— — l ^ 

+ e -(2D T+l u: )(t+r-s) Yi _ i ^ e -(l(l+l)D r + lm u: )s j (33) 

After introducing these expressions into eq. ( |3C| ) we will obtain three terms corresponding to an expansion of the 
correlation function in powers of A(i), of zeroth, first and second order respectively. The presence of this driving yields 
an explicit dependence of the correlation function on the time t, instead of depending only on the time difference, 
as occurs in the stationary case. The method for removing this dependence on the initial time is to average the 
correlation function over a period of the driving p]]. After doing this the first order term vanishes, and consequently 
we do not worry about it and compute only those whose average gives a non-zero contribution, i.e. the zeroth and 

second order terms. Taking this into account, and by applying eq. ( |l2| ) to 9(v, t\Ro, to) 



(R x (t)R x (t + T)\R (t )) ~ J dv J duu x v x y Q (v,t\R Q ,t Q )y Q (u,t + T\v,t) 

+ dv duu x v x Aty(v,t\R ,t )A^(u,t + T\v,t), 



(34) 

where the sign ~ indicates that the terms which vanish after averaging over the period of the driving has been 
neglected (although the average has not been performed yet). After introducing the corresponding expressions for 



^(v, t\Ro, to) and by using eq. (33) we obtain (the details of this computation are given in Appendix B) 

2 



(R x (t)R x (t + T)) 



|) 2 2 e - 2 ^cos(^ )+ (| X:! 



X 2 (t)e luJT 



dt'e iut Xx(t') 



(35) 
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where we have defined 



(R x (t)R x (t + r)) = — I dR (R x (t)R x (t + T)\Ro(t )). 



(36) 



At this stage, and before applying the Fourier transform to the correlation function to obtain the power spectrum of 
the process R x (t) , we average eq. (|35|) to obtain 



(R x (t)R x (t + T)) = — 

ZTT 



dt(R x (t)R x (t + r)) 



i) 



= 1 T 1 - 



2D,.t 



. 2ir 

cos(w t) + ( — 



\ 2 ILUT 



dt'e lult ' 



xS') 



(37) 



This computation has been carried out in the assumption that r is a positive quantity. To extend our computation 
to r < we have to use the backwards Fokker-Planck equation. The operator which generates the backwards evolution 
of the probability distribution — C) Jll| being C the Fokker-Planck operator and & its adjoint operator. Consequently 
the formal solution of the backwards Fokker-Planck, equivalent to eq. (E2|) is given by (t < t Q ) 



V(R,t) = e-C-'o^tf^io) + / <ft'A(t')e- (t -*' 3 ^£i*°(*0- (38) 

J t 

As in the case of the operator Co, the spherical harmonics are eigenfunctions of — C with eigenvalues given by 

- CtY lm (k) = (1(1 + l)D r - imLu )Y lm (R), (39) 

Thus the process to follow in the calculation of the correlation function for r < is identical to the corresponding 
computation for r > changing the eigenvalues of the operator Co by the ones of — Cq, which yields 



4"7T 

(R x (t)R x (t + T)) = [ — 



2e 



2D r r 



. 2tt 

cos(wqt) + ( y 



\ 2 iur 



dt'e iu}t ' 



xS') 



(40) 



for t < 0. 

We now apply to this averaged correlation function (now defined for — 00 < r < 00) the Wiener-Khinchine theorem, 
which states that the power spectrum and the correlation function are related through a Fourier transform. Thus, 

/■oc 



Q(n) = 

N(to) = 



dr (R x (t)R x (t + T))e lUT = N(to) + S(oj)5(Q - w) 
2 



2D r 



2D r 



AD 2 r + (Q + Lo ) 2 4L>2 + (O _ WQ ) 



(41) 



Since our purpose is to discern whether or not SR is present in the relaxation process of the quantity R x (t), we 
proceed to compute the signal to noise ratio, R, i.e. the ratio between the weight of the delta function in (]4l|)i and 
the noisy part of Q(tt) computed at the frequency of the driving. From eqs. (p5|) and (Ed]) we achieve 



_ S(m) 2 6 



2D r 
4D , i + (ui+uJo) 2 



2D r 
4£)j? + (w-w )' 



~'+^o 



41>2 + (w+w )2 



4D 2 +(uj-ujo) 2 



2D r 



2D r 
4Dj!+(tj-wo) 2 



(42) 



This quantity has been plotted in Fig. 1 as a function of the inverse of the Peclet number, Per 1 = D r /u)o, which 
measures the ratio between the time scales associated to diffusion (thermal noise) and flow. The presence of a 
maximum in R for a non-zero value of this parameter shows the existence of stochastic resonance in the relaxation 
process of a dipole in a shear flow. In addition to the slow relaxation to the single attractor of the dynamics, our 
model includes another effect, which hides, to some extent, the SR profile. To understand this, note that even though 
the signal is too weak, it nevertheless causes the position of the attractor of the dynamics to vary, and so the output 
will always have a nonzero component at the signal frequency. This fact causes the SNR to go to infinity in the zero 
noise limit |l9|. 
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V. MEAN FIRST PASSAGE TIME 



In this section we study the behaviour of the escape time distribution and the mean first passage time of the 
magnetic dipole immersed in a shear flow. To this end, we have to account for the fixed point orientations of eq. (Q) 
in the case Aq < 1. In this situation there is a single fixed point corresponding to an orientation contained in the 
plane x — or, equivalently, (f> = ir/2. However, when X(t) > this stationary orientation is in the subspace z > 
(cos# > 0) and in the subspace z < (cos 9 < 0) if X(t) < 0. Therefore, we are going to study the escape from the 
region z > (cos 9 > 0) assuming that the initial orientation of the dipole is contained in this region. Consequently, 
we have to solve the Fokker-Planck equation (|l0|) with absorbing boundary conditions in the plane cos 9 — |^0|,|lJ , 



*(cos6> = 0,0, i) = 0. 



(43) 



Since this escape problem will be treated perturbatively, the first step is to analyze the eigenvalue problem of the 
operator £q eq. ([ll]) under the boundary condition ( ff3| ) . It is easy to check that the eigenfunctions and the eigenvalues 
are the same with the restriction that only those spherical harmonics which vanish at cos 9 = are solution of this 
eigenvalue problem. From the parity properties of the associated Legendre functions Jl5| one can see that (43) selects 
only the spherical harmonics such that I + in = 2n + 1 with n — 0, 1, .... Thus, we have 



tf(JU) = J2 aim(t)Yi m (k), 



(44) 



(2,m) 



where (I, m) denotes that the sum is carried out over < I < oo and —l<m<l restricted by I + m = In + 1. 

In order to evaluate the mean first passage time (MFPT), we have to compute previously the survival probability, 

S(Ro,t), and the escape time distribution (ETD) which are related through 



p(Ro,t) = - 



dS(R a ,t) 
Jt ' 



(45) 



where S(Ro,t) is defined by 

S(k ,t) 



■R 



2tt 



dR^{R,t\R )= d(f> d(cose)^(cose,4>,t\R ), 



(46) 



with 1Z, the region from which we are studying the escape problem (in the present case cos 9 > 0), and Rq S 1Z, the 

initial orientation of the dipole. The probability distribution ^>(R, t\Ro) is obtained from eq. ( |l2| ) with the boundary 
conditions (43) and the initial condition 



V(R,t = 0) = 6{R-R ) 



(l,m) 



YC m {R )Y lm {R). 



(47) 



Before proceeding to obtain the survival probability, there are some facts to consider which will facilitate further 
computation. Looking at eq. ( (46| ) one can realize that, due to the integration over the azimuthal angle, only terms 

with m = contribute to S(Ro,t). Consequently, the selection rule I + m = 2n + I reduces to keep only the odd 
values of I. In addition, we are interested only in the modes with greater relaxation times. Therefore, from the whole 
series eq. ( pf4| ) we are only interested in the term I = 1, m = 0. Thus, our purpose is to obtain the coefficient aio(t) 
up to first order in \(t) from eqs. (10) and (|43|). Up to zeroth order, we have 



.(°) 



(*) 



-2D r t 



Y W (R 



(48) 



The obtention of the first order contribution, a^g (t) , requires a little bit more elaborated calculation. To proceed 
further with this computation, the operator C\ acting on Yi m (R) yields 



&Y lm (R) 



-2uo\l y (lu (#) + Y 1 - 1 {R))Y lm (R) 



w 



—Y w {R)TZ y - iJ — (Y n (R) - Yi^ 1 (R))TZ 2 



Yl m (R), 



(49) 
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where the action of 1Z Z on Yi m (R) is given by eq. ([l4"l)i and 



n y Yi m (R) = -\{V(l - m)(l + m + V)Y lm+x {R) - ^(l + m)(l-m+ l)Y lm -i (R) } ■ 



(50) 



From eqs. ([49]) and ( p0|) together with the rules for the addition of angular momenta familiar from quantum mechanics 
[ filj and the selection rule / + m = 2n + 1 imposed by the boundary condition ( |4S| ) , one can deduce that only the 

terms Y 2 ±i(R) contributes to a^ (t). The mentioned rules of addition of angular momenta implies that the product 
of two spherical harmonics, Yi m Y pq has projection onto a third spherical harmonic, Y rs , only when m + q = s. On the 
other hand these same rules imposes the restriction that the product Yi m Y pq only projects onto the subspaces such 
that \l — p\ < r < I + p. By using these restrictions one can see that when in eq. ([49J) one takes I = 1 and m = one 
obtains a vanishing contribution and only when I = 2 and m = ±1 the contribution to (t) is different from zero. 
All other contribution of higher values I are explicitly excluded by the rule \l — p\ < r < I + p. 
Taking these considerations into account and by using the following results 



47T 



4tt 



4tt 



2n 



d<j> / d(cose)Y 11 (R)Y 2 - 1 (R)Y w (R) 
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80^- V 3tt 



In 



^ d(cose)Y 1 ^(R)Y 21 (R)Y 10 (R) = ^, 
i „ 



# / d(cos6)Y 10 (R)Y 20 (R)Y 1Q {R) 



3 / 5 
20tt V 4?r' 



(51) 



we obtain the first order correction to the coefficient 



a< >(t) 



aw 



(t). 



Ap [15 f V2 A \ g 4.P r + z(w - trip) , v ,g , 4P r + t(tri + irio) 
^Q-\/ — ( — +2V3 \\ Y 21 (R ) ^ + y 2 _ 1 (ie )______ 



40 V 4tt 



16£)2 + ( w - Wo ) 2 



16£>2 + ( w + wo) 5 



(52) 



Eqs. 



and (p2J) together with (kt4J) and Mfl) allows us to obtain the survival probability S(Rq, t), which is given by 



S(R Q ,t) = n x — Y w (R )e 

47T 



-2£>,t 



+y 2 _i(i? )e~ 2I?rt / dT\(T)e-( iD "-^ T 
Jo 

This quantity is directly related to the MFPT, since 
where we have used eq. (ffq). Consequently, the MFPT is given by 



drA(r)e- (4A ' + ^ o)r 



dttp(R ,t) = / dtS(R ,t), 
(l Jo 



(53) 



(54) 



T(i? ) = T (R ) + AT(R a 



T (Ro) = 7T- 
AT(i? ) = T 



3 yT^A2 
4^ 2D r ' 

40 V 4tt I 2 / °V 24tt \l6D? + (uj - uj ) 2 WD? + (cu + iv ) 2 



(55) 
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where we have taken i? — R s - I n Fig- 2 we have plotted the quantity AT/T a . This figure shows that this quantity 
exhibits a minimum, as required for the appearance of SR. 

The knowledge of the survival probability allows us to obtain the ETD, p(Rq, t). From cqs. ( f45| ) and ( |53"| ) the ETD 
is given by 



p(R ,t) = e 



-2D r t 



Y 2 -i(R ) 



2D r iT\ —Y w (Rq) - A ( 

V 47T 

2D r + + luq) 
16L> 2 + {uj + lj ) 2 



Y 2 -i(Ro)e 



-(4D r +i(uj+u ))t 



Ac, 



2.D r + i(u> — luq) 



_L_ /l5 / V2 
40V 4^ i T 



16L> 2 + ( w _ Wo )2 

2V3 ) {y 21 (^ )e-( 4D " +l ^-" » t 



(56) 



By taking R = R s , we finally obtain 

40tt 



-4D r t 



16D2 + - Wo )2 16Z? 2 + (w + w ) 2 



-(sin[(o; — wo)t] + sin[(tj + u>o)t\) 



(57) 



where Ap = p — po, being po the corresponding ETD when the amplitude of the oscillating field is set to zero. 

The succession of maxima in the ETD (see Fig. 3) is indicating that the dynamic of the magnetic dipole suspended 
in a shear flow under a periodic field exhibits SR. 



VI. DISCUSSION 



We have shown that the relaxation process of a dipole immersed in a shear flow exhibits SR upon application of 
a weak periodic field. To this end we have computed three quantities typically used to characterize SR, namely the 
signal-to-noise ratio, the escape-time distribution and the mean first passage time. All of them behave as expected 
for a process in which SR occurs. 

Previous works devoted to analyze whether or not SR is present in the relaxation process of an overdamped dipole 
in a fluid at rest have revealed that this phenomenon does not occur in the linear regime JT||]. Effectively, linear 
response theory (LRT) predicts a maximum in the signal, i.e. in the susceptibility, as a function of the noise level. 
Unlikely, the SNR decreases monotonically with the noise level. This behaviour can be easily understood. In the limit 
of zero noise the output of the system has a small component (proportional to the applied field) at the frequency 
of the signal whereas the background noise vanishes at zero noise level, being this behaviour the responsible for the 
monotonic dependence of the SNR on the noise intensity. 

In our case, the situation is completely different. When the fluid in which the dipole is suspended is submitted 
to a pure rotation (vortex flow), both output signal and output background noise exhibit a peak at the same value 
of D r (see Fig. 4). Consequently, although the background noise vanishes when D r goes to zero, the characteristic 
SR profile of the SNR can not be completely hidden, as shown in Fig. 1. This feature arises as a consequence of 
the presence of shear acting on the suspension, thus, the appearance of SR in the system studied in this paper is a 
non-equilibrium feature. 

In a sense, the mechanism yielding SR in this system is similar to the one operating in SR in threshold devices ||,[l0) • 
Due to the presence of noise, the dipole can eventually acquire enough energy to get out from its stable orientation 
by crossing the absorbing barrier (the threshold) cos 9 = 0. After this, the system is driven to its stable position. 
This process produces a short spike in the magnetization. Of course, the time that the system takes to return to the 
fixed point has to be smaller than the semiperiod of the oscillating magnetic field. Thus, SR in this system can be 
understood in the same way as, for example, the SR in level crossing detectors H. 

LRT has been one of the widest used tools in the study of stochastic resonance E2] . When the system, in the absence 
of the external periodic force, is in thermal equilibrium a very adequate form of describing stochastic resonance is in 
terms of the susceptibility. This is because the noisy part of the power spectrum is given directly by the susceptibility 
through the fluctuation-dissipation theorem, 

lm A(^) = tStN(Q), (58) 







This result is correct when the fluctuations whose spectral density is given by N(fl) have the thermal equilibrium 
state as reference state [H. 

However, in the present case we are dealing with a system which is maintained in an out-of-equilibrium steady 
state due to the presence of a shear flow. It is evident from Fig. 5, where we have plotted the imaginary part of the 
susceptibility corresponding to R x and the noisy part of the power spectrum, that these two quantities are clearly 
different. Note, that if loq = 0, i.e. the system in the absence of the periodic field is in equilibrium, the relation ( [58|) is 
fulfilled. Thus, we have shown that although we can define a susceptibility which describes the response of our system 
to a small perturbation we can not describe SR by means of LRT. The reason can be found in the fact that due to 
the non-equilibrium nature of the attractor of the dynamics the fluctuation-dissipation theorem fails to be valid. 
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APPENDIX A: LINEAR STABILITY ANALYSIS OF THE FIXED POINTS OF EQ. (§) 

From eqs. (|J) and (f|) one can see that the time derivative of the R vanishes either when f2 p = or when £} p x R = 0. 
In the first case we have 

(l p — lo jz + A (^R s X xj I = => z = -X R S x x. (A.l) 
From this equation, and taking into account that \R S \ — 1, we obtain that the stationary orientation is 



R s =±^l-X 2 x + X^y. (A.2) 

This solution exists only when Ao > 1 and corresponds to a fixed orientation of the dipoles, given that the intensity 
of the magnetic field is high enough to maintain this fixed direction. 
The second possibility leads to 

zx R s = -X a (R s x x) x R s = -X (x- (x- R S )R S )- (A. 3) 



Eq. (A. 3) provides two equations for three unknowns. If one sets R z = one recovers eq. (A.2). If, by contrast one 



makes R x = then a new stationary orientation is obtained, 



R s = X y± yjl - X 2 z, (A A) 

which exists only when Aq < 1. This orientation gives rise to a rotation of the dipoles with angular velocity 



n s = u ^l-\ 2 ^l-\tz±\ yj, (A.5) 

since, in this case, the field is not strong enough to inhibit the rotation caused by the shear flow. 

The linear stability of this fixed points is better analyzed in spherical coordinates. Taking into account that 

(R x ~x) x R = S- R{R- £), (A.6) 

we obtain 

•M 1 - r2 x) - Ry> 
— XR X R y + R x , 

-XR X R Z . (A.7) 



1 


dR x 


w 


dt 


1 


dRy 


Uq 


dt 


1 


dR z 



wo dt 
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After expressing the components of R in spherical coordinates, we obtain the following bidimensional dynamical 
system 

— = A cos cos <p, 

too at 

-§ = -A^ + l- (A.8) 

luq at smO 



where 9 and <j> are the polar and azimuthal angles, respectively. By linearization of eqs. flA.q ) around the A > 1 fixed 
points we obtain the following matrix 



A(A>1)= V " XJ , — T (A.9) 




which implies that, if A is positive, the orientation corresponding to choose the sign "+" in eq. ( A. 2 ) is stable, while 
the other one is unstable. 

The same linearization procedure carried out around the A < 1 fixed points leads to 



A(A<1)=( , 2 .2-^ AV !r^) <A.H» 



-\ 2 Vl^ A2 
The eigenvalues of this matrix are given by 

a = ±iAVl-A 2 . (A.ll) 



APPENDIX B: COMPUTATION OF EQS. (|33|) AND (g^). 



In this appendix we work in detail some steps of the computation of the power spectrum corresponding to the 
relaxation process of a dipole under an oscillating magnetic field in a shear flow, in particular we calculate the 
integrals which yield eqs. (33) and (p5|). From eqs. ( |3l| ) and (|32|) 

duu x 'fy{u, t + t\v, t) = Ii + ^2) 

duu x ^o(u, t + t\v, t) — j duu x e TC °S(u — v), 

t + T 



h 



h = J dUu x I^{U,t + T\%,t) = J dUu x J dsX(s)e it+T - s)Co £ 1 e sCo S(il- : u) 



To begin with we focus on the integral 7i, which can be rewritten as 



du ( e T »u x ) 8{u — v), 



where £q is the adjoint operator of Co defined by 

J ^A(£ B) = J dtt(4A)B, 
with A and B two arbitrary observables. Explicitly, Cq is given by 

C\ = uj R -TZ + D r K 2 , 

C\Yi m {R) = (-1(1 + l)D r + ioj m)Y lm (h), 



(B.l) 



(B.2) 



(B.3) 



(B.4) 
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and therefore eq. (B.2) reads 



/X = -Jy / ,/,7 



e -(!U3 r -i U o)Ty ll ^ + e -(2fr+^o)r Fi _ i ^ 



leading to eq. (J33|)i. In eq. (B.5) we have used the relation 

f2w 



To compute the integral J 2 ? w e have to use the following representation of the delta function 

oo I 



1=0 m=-l 



After introducing this expression into eq. (B.l)3 we obtain 

oo I rt-\-T 



1=0 m=-l Jt J 

°° 1 r t+T r / ♦ \ 

= E E Y £»W J ds\( S )e-W + V D ^^ j du^^u^CMn), 



1=0 m=-i 



-(2Dr-iu )(t+T-s)Y 11 (g\ 



and by using eq. (B.6), eq. (B.E) yields eq. (|33|), ie. 

/ 2 = -i/ T E E Y L@) / / 

1=0 m=-l Jt J 

Once these expressions have been obtained we can compute the correlation function given by eq. (|3i 
{R x {t)R x {t + T)\ko{t )) ~ / d0» 1B *o(fl r ,t|^o,*o)i'i+ / d«>xA*(lU|^o,*o)l2, 



being # (u,i|-Ro,*o) and A#(tf,i|.Ro,*o) given by 

M*A$o) = e tC "5(t-tio) -EE ^(^o)e~ (/(i+1)A ' +m ^ o)t l, m ^), 



1=0 m=-l 

dr\{r)e {t - r)Co C 1 e rCa 5{ij'- R ) 



o 

oo I 



1=0 m=-l 



E Y i*n(Ro) / drA(r) e -w i + 1 )^+™- c ') r e ( t -'-) c «/: 1 r im (^). 



(B.5) 



(B.6) 



(B.7) 



(B- 



(B.9) 



(B.10) 



(B.H) 



where t he init ial time to has been fixed to zero and eq. (B.7) has been used. 

Eqs. ( B.ll ) provide the evolution of the probability distribution under the condition of the system being initially 

in the state i?o- Since a priori nothing is known about this initial condition, we assume that R is a random variable 
uniformly distributed over the orientation space, consequently we average the correlation function over the distribution 
of initial states (eq. (|36|)). Taking into account that 
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— / dRoY* m {Ro) = SiflSmfli 

the average of the correlation function over initial conditions 

{R x {t)R x {t + t)) ~/ 3 +/ 4 
h = J dvv x Ii, 

dvv x f dr \{r)e^ Ca dYo^h. 



After introducing eqs. (B.5) and (|B.q) into (B.13)2, we obtain 



h = 



1 

4-7T 



d&y (HifS) + n-i(£)) [e"^-^)^!^) + e-f^+^^y^x^) 



-2D r r 



cos(ujqt), 



where the orthogonality relation for the spherical harmonics, 



rtv*rsr\v rZr\ 4tt ([+™)! 



has been used. On the other hand from eqs. (B.9) and ( B.13| ) 
1 /2tt 



VT^ / drX ^ / dsA ^ / ^{(^ii(^) + n-i(^))(e- (2A '-^ )( *- r) rii(^) 

+ e-( 22> ' +<wo ^ t - r )yi-i(l0) / ^{[e-( 2 ^-^)( t+T - s )Fii(S) + e-^+^)(* +T - s )yi_i(5) 

e -a(i+i)^+im W0 ) S £ lFjm ^Jl 

Let us focus our attention on the integral over v. 

d0(Yn@) + Fi-i(tf)) (e-< 2D -- iwo >(*- r 5yii(l0 + e- (2 ^ + ^ )(t - r) ri-i(^)) ^(5) = 

^^(^e-^-^^^^yn^y^^) + / d^n^e-^+^-^^Y^). 



(B.12) 



(B.13) 



(B.14) 



(B.15) 



(B.16) 



(B.17) 

From the rules of addition of angular momentum familiar from quantum mechanics, which imply that the product 
of two spherical harmonics, Y pq Y rs , only has a non- vanishing projection over a third spherical harmonic, Y, m , when the 
relations \p— r\ < I < r + p and q + s = m are fulfilled, it is easy to see that these integrals will give a non-zero result 
only when I = 0, 1, 2 [Q. In addition, for the integrals containing the products Yi±i(v)Yi±i(v) the parameter m has to 
be m — ±2 whereas it must be m = for the integrals with Yi±i(v)Yi T i(v) to yield a non-zero contribution. However, 
although these integrals gives in principle a non-vanishing contribution, note that when we perform the integral over 
the variable u in eq. ( B.16| ) the terms introduced by these contributions finally yield, by the orthogonality property 
of the spherical harmo nics, a vanishing result. Thus, only when I = and m = contributes to / 3 . Taking this into 
account and using eq. ( B.15 ) 



3 r t 



dr 



t+T 



ds 



A(r) | e -( 2D r-iwo)(*-r) + e -(2D r +iu> )(t-r) j 
A(s) | e -( 2D r-iw )(t+r-s) + e -(2D r +iw )(t+T-s)| 



(B.18) 
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finally performing the changes of variables t' — t — r and t" = t + r — s and using eq. ( |2q ) we obtain 
In this integral the upper limit goes to infinity by causality reasons 



h = ( ¥ ) A a (i)e- toT ( / dte-^ Xx (t) I . (B.19) 
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FIG. B.l. Signal to noise ratio as a function of the inverse of the Peclet number Pe 1 = D r /u>o- We have represented 
non-dimensional quantities. 
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FIG. B.2. Mean first passage time as a function of the parameter a = u>/u>o- The presence of a minimum reveals the existence 
of stochastic resonance. Pe~ = D r /u)?> = 0.08. We have represented non-dimensional quantities. 
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FIG. B.3. Escape time distribution for a dipole in a shear flow under an oscillating magnetic field. The succession of maxima 
is a signature of the presence of stochastic resonance in this system. Pe _1 = D t /ojq = 0.08 and a = w/wo = 0.7. We have 
represented non-dimensional quantities. 




FIG. B.4. Output signal and output background signal as function of the inverse of the Peclet number. Solid line repre- 
sents the quantity iS(.Pe _1 )u;o/Ao whereas dashed line represents iV(Pe -1 ). We have taken a = 0.1. We have represented 
non-dimensional quantities. 
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FIG. B.5. Comparison between the imaginary part of the susceptibility of the signal R x (solid line) and the noisy part 
of the spectrum computed from the Fokker-Planck equation (dashed line). Pe^ 1 = D r /u)o = 0.08. We have represented 
non-dimensional quantities. 
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